State  University  of  New  York  at  Buffalo 
Department  of  Computer  Science 


Research  Office.  (('rani.  DA  AC> 


Unc  la  s sifted 

SLCUMI  ' Y Cl  ASSlI  ICATION  Of  THIS  I’Al.U  fnl.rtdj 


REPORT  DOCUMENTATION  PAGE  KorkDcompu  Tnlcroim 

i lioMiu.K  ” oovF  AC  Cl  SSIOW  NO-1  j.  KCClflENT't  CATALOG  NUMtttH 


T | T t l (*nd  Oiibtltl/') 


Techniques  of  Quantile  Regression 


1.  All  1 HOHi'lJ 


S lyri  Of  HI  I'OHl  A I LHlOO  COVtlNt  O 

Technical 


. HCRf  ORMINC  01(0.  R!  PORT  NUMBER 


II.  CONTRACT  OR  CHANT  NUMB  R(«J 


Jean-  Pierre  Carmichael 


DA  AC.29-76  Cl-0239 


nil  ORMINO  OHO  ANI  2 AVION  N AMI  ANOAUORfSS 

Statistical  Science  Division 
State  University  of  New  York  at 
Amherst.  New  York  14236  

It  CON  I ItOLUNC  Cl  f I l<  f lltMl  ANDyAOONESS 

II.  S Armv  It'  e irrh  Oilice 
l-'or.l  Ollier  hot  12211 
Research  I'tutk’lo  I'.iilc  NC  27/09 


111.  I'nOOCAM  M I Ml  tJ.I'ROU  T.  TASK 
AMI  A ft  KiftK  UN  n NUM'.ltPi- 


\:  R ( PORT  D AT  b 


June  1978 


U NUMiti  It  C»  r AGKS 

27 


\ MON  I RING  ACtNCV  NAME  A Al  Rl  II  1 ■ / . i Controlling  Ottleo)  15  SECURITY  CUAS&.  (o<  IM« 

tlm  las  si  ried 


l;-«.  I ! n A ••  I ' |C  ATION  OOMrNCHADINl. 
SCHtOULl 


lo.  01  SO  III  HU  llON  !.T  AT  TMI  N 1 («l  Ihl . Kcpnrt) 


Approved  for  public  release;  di  r»t  ribution  unlimited. 


17.  DISTRIBUTION  f-1  Ail  Mi  NT  (■>!  r ’■  mb  ptrmrt  It:  /»Mclr  20,  H dlllmrmnl  from  Krpotl) 


18  51'PPl.LMl  N1  AMY  NOUS 

The  findings  in  this  report  are  not  to  be  construed  as  .ht  official 
Department,  of  the  Army  position,  unless  so  designated  by  other  authorized 
documents.  


1£>.  MY  V OhP'i  (Continue  on  r<*voi.-o  ml  Jo  II  n onrl  Identity  by  bU>  <k  numb*-) 

Density  Estimation 

Nonparametric  Regression  % 

Order  Statistics 


20  ABSTRACT  (ContloTirr^tg^f^o  mldm  II  n memomory  mud  I'mntlty  by  block  nunbrnr) 

We  study  Jhe  asymptotic,  properties  of  different  nonparametric 
estimators  of  a regression  function^.  The  motivation  for  these  estimators 
comes  from  the  unified  approach  to  statistics  developed  by  Parv.cn  under 
the  name  of  ’"flonpo  ramet  r ic  Statistical  Data  ScicMceM'Tn  which  the  quan- 
tile function  plays  a crucial  role.  We  call  them  quantile  regression 
estimators,  ,— <• 
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TECHNIQUES  OF  QUANTILE  REGRESSION* 

by 

Joan- Pierre  Carmichael 


Intrude c t i on 

Given  observations  {(X.,Y.),  i=  l,...,n)  on  random  variables 

(X,  Y)  with  joint  distribution  F v(x,  y)  , we  want  to  estimate  the 

X,  i 

regression  function  of  Y on  X , E[Y|X  = x]  , nonparametrically. 

In  order  to  find  a natural  estimator  (simple  computationally  and 
intuitively  appealing),  Parzcn  (1977)  developed  the  following  theoretical 
approach. 

1.  Theoretical  Approach; 

Let  U = F (X)  and  U - F (Y)  , then  the  joint  distribution 
1 A ^-1 

of  Uj  and  U^  is 

E>Ul,U2(urU2)  = lX,  y(QX(V’  °Y(U2)) 

and  their  joint  density  is 

fx.  y(°X(u1  *’ 

U.'U2  ‘'“2  fx(QX<Ml»  hVl') 

where  F is  the  distribution  function  of  Z 

/a 

f is  its  density  function 

fit 

Q^,  is  its  quantile  function 
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Let  r(x)  be  the  regression  function  of  Y on  X - x . 


«,  y fx  Y(x,y)  dy 


r(x)  = E[Y|X  = x]  - J 


We  now  define  the  regression-quantile  function  rQ  (•)  by 


•Q(u)  r(Qx(«))  ==  E[Y|X  = Qx(u)] 


llow  do  we  compute  rQ(»)  ? 


By  definition. 


. „ r yfx.  Y(°x(ul-ydy) 
fx(ux(u)J 


Let  y = Qy(u?)  , then 


rQ(u)  = JQ  Qy(u2)  dy  tU  (u.  «-,) 


r 2 


If  wc  introduce  a Dirac  delta  function,  we  can  express  rQ(*)  as  a 


double  integral 


rQ<u>  = It!  Io  °Y(V  6(U1  • u)  d DU1,U2(U1,U2) 


We  estimate  rQ(»)  by 


* Jo  fo  W h(^  K (\^-)d  %,.U<U1-V 


r 2 


3 


function  of 
function. 


U j and  U.^  . 


is  an  estimator  of  the  joint  distribution 
It  could  be  the  empirical  joint  distribution 


cyo.  is  an  estimator  of  the  quantile  function  of  Y . It  could 
be  the  empirical  quantile  function  of  the  Y's  . 

K(*)  is  an  approximator  to  the  Dirac  del'.a  function. 


2.  Different  Estimators: 

Let  > be.  the  observation  associated  with  X|^  , where 

X.,.  < X,~,  <•  . .<  X.  . . Yr.  i is  called  the  concomitant  of  the  i*^ 
(1)  (2)  (n)  U : n] 


order  statistic. 

A 

Let  D (•  , •)  be  the  empirical  joint  distribution  function 

uru2 

of  Uj  and  . It  has  jumps  of  si/e  1/n  at  points  of  the  form  (i/n,  R. /n) 

where  R.  is  the  rank  of  Yr.  . among  the  Y's  . 
i [i  : n] 

Let  I<  (•)  be  a -kernel  function  with  bandwidth  parameter  h(n)  . 

A 

A first  estimator  of  Q^(*)  is  given  by 


Q1<V  = Y[i  , n] 


s 


u„ 


Then,  equation  1.2  becomes 


2.  1 


a n 

rQj(u)  = 


j/n 

Y[j:n] 


n 


h(n) 


K 


dt 


4 


Usually,  as  in  Yang  (1977),  this  is  approximated  by 


2.  2 


rQ 


Yang  studied  the  statistical  properties  of  that  estimator. 

A 

Note  that  rQ^^(*)  can  be  viewed  as  the  result  of  smoothing 

amplitudes  Yjj  t observed  at  equidistant  points  of  the  form  j/n  . 

Clark  (1977)  recommends  to  interpolate  linearly  between  the  successive 

points  {(j/n,  Yr,  .)}  to  get  an  estimator  with  maybe  more  derivatives 
U ! “] 

than  the  kernel. 


Define 


W 


[Isn]  * 


0 £ u_  s 1 /n 


Y(  j : n]  ('  + 1 'nU2)  + Y[j  + l :n]("U2-j)’ 


J * u_  s 
n 2 


j_+  1 


j = 1 ..... n - 1 


Then 


2.  3 


•V'>  = Jo 


It  lias  been  remarked  before  that  it  is  difficult  to  smooth  a curve 
unless  it  is  relatively  flat.  People  would  then  recommend  to  subtract  a 
trend  term  from  the  data  before  smoothing. 

We  would  like  to  propose  instead  to  smooth  the  first  differences 


Q#2<*>  * 


...  --SOS 


5 


°2  (U2}  = 


0 s u2  < 1/n 


" ' (Y[jtl:n]'Y[j 


< ~ J + l 

S«,  <* 

2 n 


j = 1 , . . . , n- 1 


We  then  form  tin;  estimator  rQ  (•) 


*ciw  ■ fo1  “'2“>  1^) 


rOjIu) 


-ni  ,-Yt.  ^rh  fj/1/nKfT7-?U 

•_1  \ [j+ 1 : n]  [j  : n]/ h(n)  Jj/n  \ h(n)  / 


A A A 

Z.S  r03(u)  = JjU/2  rcy  (s ) ds  4 rO^l/2) 


Because  an  estimator  of  rQ(*)  would  be  the  indefinite  integral  of 

A A A 

rQ'  (•)  , we  fix  the  value  of  rQ^fl/2)  to  be  rQj(l/2)  as  we  feel  that  all 
estimators  are  usually  good  for  the  middle  values.  The  problems  and  the 
differences  between  estimators  usually  appear  near  the  endpoints. 

A 

Finally,  we  can  smooth  Q (*)  using  the  autoregressive  method 
by  computing  its  Fourier  coefficients 


<P  (v)  = e2TTltV  Q 2(t)  dt 


j+1  /n 


cn(v)  = E n(yr  - Yr . A J 


2rt  itv 


| v | = 0,  1 , 2,  . . . , m 


f 
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From  the  CP  ( • ) ' s , we  compute  the  autoregressive  coefficients 
by  solving  the  Yule- Walker  equations 

A 

2.7  rQ^(u)  = 

and 

A A A 

2.8  rQ4  = JjU/2  rQ2  (s)  ds  + rQ^l/2)  . 

Note  the  relation  between  the  linearized  version  of  the  data  and 
first  differences.  Taking  k order  differences  would  be  like  interpo- 
lating between  data  points  with  a k**1  degree  polynomial. 

3.  Statistical  Properties: 

3.  1 General  Results 

Yang  (1977)  studied  statistical  properties  of  linear  functions  of 
concomitant  of  order  statistics. 

Among  the  different  estimators  proposed  in  the  previous  section, 

A A 

only  rQ2(*)  and  rQ^(*)  are  not  of  that  form. 

For  convenience,  we  reexpress  the  three  major  results  of  Yang  in 
a form  lnore  related  to  our  purpose. 

We  need  the  following: 


k „ 


1 + 


a. 

j. 


2tt  ij vi  | 2 


7 


M„  ■ Jo1  Jo1  «<v  iii K 


iii  K(n^r)da.1,u2(ui'u2l 


where 


g(u  ) = n(x...,  Yf.  , AjlJ.  *u,  <- 
2 \ (l)  [ i : n]  / n 2 n 


CX(x)  = E[1I(X,  Y)  | X = x] 


o (x)  = Var  (I1(X, Y ) X x) 


Assumptions 


A1  - E[  |11(X,  Y)|  J < » 

A2  - q(x)  can  be  expressed  as  a difference  of  two  increasing 

right- continuous  functions 
l 

A3  - 0 (x)  has  the  same  property  as  ci(x)  or  F,(x)  is 

A 


absolutely  continuous 


<x(Q(t))  i 


is  continuous  at  t - u 


A5  - E[H(X,  Y)  1 < « 

A6  - a'  (x)  = ■—  Ci(x)  exists  and  a.'  i^Q(t)y  is  continuous  at 


t - y , 0 < u < 1 

A 7 - — — a 

dt 


M 


exists  and  is  continuous  at  t = u 


131  - There  exists  M > 0 such  that 

I K(tj ) - I<(t2)1  < M • |tj  - t?J  for  all  t ^ t2 

132  - | tK (t ) I -*  0 ns  1 1 1 -»  1 
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B3  - J_j  K (t ) dt  = 1 

B4  - lim  h(n)  = 0 
n-»» 

B5  - lim  h-|(n)(1°tn1°Bn)1/4  =0 
rr» » 

B6  - J_\  t K(t ) dt  - 0 
B7  - t2  ( K(t)  \ dt  < 00 

B8  - K"(t)  exists  and  satisfies  B1  and  B2 

Thmi  - Consistency  (Yang's  Theorem  5) 

Under  assumptions  A1  - A4  and  B1  - B5, 


lim  E[MJ  = afQ(u)J 
n -» « 


9 


Th  2 - Asymptotic  normality  (Yang's  Theorem  6) 
Under  assumptions  Al  - A6  and  131  - 135  , 


^/nh(n)  (M^  - E[ -»  N ^0,  CJ2Jq(u)^  R2(t)dt^ 


Th  3 - Asymptotic  bias  (Yang's  Corollary  1 to  Theorem  6) 
Under  assumptions  Al  - A7  and  13 1 - 138  , 

E[m1-Q;W  AZ  , x 1 , 


1 n*  "S, 
lim 

n-*®  h (n) 


- - a(o(u)N  • f t?K(t) 

, b \ / V 


/ , AD  ! v 1 \ 

/nh(n)  M - a Q(u)  ->  N .0,  o\Q(u);  • f K (t)  dt’ 

\U  K 'I  \ K ‘ -1  / 


Let  vis  apply  these  general  results  to  the  different  estimators  we 


presented  in  the  previous  section. 


3.2  Statistical  Properties  of  rQ^p)  and  rO  ^(*1 

A 

rQ^(*)  is  the  estimator  proposed  and  studied  explicitly  by  Yang 

A 

as  an  estimator  of  E[Y|X  = Q(*)]  . Our  rQ^p)  has  exactly  the  same 
properties  as  can  be  seen  from  the  fact  that 


/♦  y 

i_  ./♦  - u\  i f j ' u 

(n)  \ h(nVl  nh(n)  \ h(n)  , 


i - 1 * 

for  *t.  £ j/n 

n j 
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Thus,  rQj  (u)  is  a consistent  estimator  of  rQ(u)  = E[Y|X  = Q(u)]  , under 
the  conditions  of  Theorem  1,  at  the  points  of  continuity  of  rQ(*)  . 

Under  the  conditions  of  Theorem  3,  the  asymptotic  bias  is  pro- 
portional to  the  second  derivative  of  rQ(») 


For  the  kernel  we  have  been  using 


K(z)  = 


15  1 1 Z\2 

76  <*  - z > 


Z s 1 


Z > 1 


the  asymptotic  bias  is  1/7  • rQ"(u)  and  the  variance  of  the  asymptotic 
distribution  is  5/7  • Var  ^Y  | X = Q(u)y  . 

It  is  possible  to  estimate  Var  ^Y|X  - Q(u)^  by  the  same  method. 


?(v|x  = Q(u>)  - i E(yd  ;n]  - rC^u,}2^  leads' 


3.  3 Statistical  Properties  of  rQ^(u) 


We  rewrite  rQ^f*)  as  f°U°ws: 


rC^(u)  = Ij  (u)  + I2(u)  , where 


n j/n 


Vu)  = j=i  hi  Y[j:nl  * h(n)  K^h(n)'dt 

n 

n j/n 

I,(u)  = S f n • (y,.  , , - Yr.  . V - t)  • 

2 j=2  \LJ-lsn]  Id  : nJ  / n 


_L_  yjt  ~ u' 

h(n)  \ h(n) . 


11 


Note  that  I^(u)  is  just  rQ^(u)  . On  the  other  hand. 


Et(Ytj-l:n]-YD.n])I  = £ rQ(s>  ^ - s)"^1] 


and  by  expanding  in  Taylor  series 


where 


j/n  . . . 

(n  ■ t}  ‘ h(n)K(Tfr)} 


dt  = 


n 


1 


2nh(n) 


K 


R.  I * 


(Ir- A ' 

\ h(n)  / 2.  2 

' n h (n) 

k:  '“"fef)1 


1 K'{^r)- 1 K'  fe 


+ R. 
Jn 


jn  6nhj(n) 


~ < t.  < , j = 2,...,n 

n j n J 


The  bias  properties  of  rQ_(*)  are  the  same  as  those  of  rQ  (•) 

^ 1 

provided  I^*)  contributes  only  to  high  order  terms. 


E[I2(u)]  = Jj(u)  + J2(u)  + R 


We  look  only  at  Jj(») 


Jj(u)  = £ rQ(s)  Zjz nh(n))_1  ) d[(py)  sj_1(l  - s)n_j+1] 


j=2 


-12- 


I 


By  Bernstein  approximati6n. 


i 


For  0 < u < 1 , 


1 


I 


rQ(s)  h (n)  K 


[<'  ( 

V h(n)  J 


ds 


u+h(n) 

= J rQ(s)h 
u-h(n) 


ds 


because  K ( • ) is  defined  only  on  (-1,1)  and  upon  integrating  by 
parts,  this  is 

h(n)  • J rQ'(u  + th(n))  K (t ) dt  = h(n)  rQ'(u) 


On  the  other  hand. 


£ l^1  rQ(s)  d(l  - s)n|  < | Jq  rQ(s)  d ( 1 - s)  | = A 


i IJq1  rQ(s)  d en|  < l/J  rQ(s)  ds  | = B 


and  for  h(n)  < min  (u,  1 - u) 


A 

h(n)* 


f 


-13- 


So, 


rQ 


goes  to  zero  as  n 

Thus,  for  0 < u < 1 


JZ(' 


(jut's  to  r.e ro  faster. 


lim  E[rQ_(u)  - rQ(u)] 

h2(n) 


lim  E[  rQ  (u)  - rQ(u)] 
n-+«°  1 

h‘-(n) 


At  the  endpoints,  the  limit  duos  not  exist. 


3.4  Statistical  Properties  of  r Q ^ ( • ) 

A i 

We  start  by  studying  rQ^  (•) 


rQj(u) 


\ 

"1  ' 


h(n) 


ri/n 

J k 

j-l/n 


<feS) 


dt 


lty  Taylor  series  expansion. 


r 


IVrV...  ]]■-{*«•>  4.!.y,‘*--'-j+1]- 


IrQj'w  = 

i—  — 


-J*  rQ(s)  2 s h(n)~  1 k' 


1 „ n 


dii  1 

u 


+ —±-  K'  — 


h(n)  / ' , 2,  v 

' ' ' 2nh  (n) 


fe)  «j-‘  <«  - 


By  Bernstein  approximation, 


fVu>l  * -;rQ(3)h-2(n)K'(i^)dS  ‘ J «')  —J-  K"(l^)ds  + B 

L J 0 0 2nh  (n) 


For  0 < u < 1 , 


; rQJ ' (u)  = -h-1(n)  rQ^u  + th(n)^K(t)  + J rQ  (u  + th(n))  K(t)  dt  + O (n"  1 ) 

. J “ 1 • 1 


Thus 


fA  ■ 1 

:[rQ,  (u)j  = 


rQ  (u) 


E rQ  (u)  - rQ  (u) 


J t2  ■ K(t)  dt  . 


4.  Case  of  X fixed 

We  study  only  the  case  where  the  x's  arc  fixed  and  equidistant 
on  the  unit  interval,  of  the  form  [j/n].”0  . The  model  is  of  the  form 

Y = f(x)  + G , where  the  e's  are  uncorrelated  errors  with  mean  zero 
and  constant  variance. 

We  limit  ourselves  to  only  two  estimators: 


<“>  ■ ski  Li  Y<tM  • K(ik>) + = Y(i/n’  • K(: 

1 • • K(w) 


/n  - u' 


and  the  estimator  based  on  first  differences 


= / f,  (■»)  ds  + f ,o/2) 
6 1/2  1 


~ I 

where  fj  (s)  = S 


^)-  V(i/n)]  • ^ k(^) 


and  Y(j/n)  is  observed  at  x = j/n 


4.1  Statistical  Properties  of  i^  ( • ) 


(v“H  • [i  f(0Kki)  +kt,j/n)  Kkikr“)  * i tlM  • KC0 


4.2  Statistical  Properties  of  f2(») 


P 1 Inv'.‘  ffe1)'  . 

L1  J = n j = o 1/n  * h(n) 


l 


Kfe)- 


f (s) 


n-*“ 


The  asymptotic  bias  is  computed  as  in  section  3 


f,  (s)  - f (s) 


h2(n) 


2 

f"'(s)  • J t K(t)  dt 

n -*  00  2 - 1 


Thus, 


u 


E|f_ (u)J— >J  f (s)  ds  + f(l/2)  = f(u) 


1/2 


and 


f2(u)  - f(u) 


h2(n) 


+ f"(l  /2)J  • J*  t2  K(t) 
^ / -1 


dt 


One  can  write  an  exact  expression  for  the  variance  of  f^*) 


Var^J*  f^sjds  + ^(1/2^  = Var^f  ^ ' (s)  dsj  + Var^f J (1  /2)j 


/,u 

+ 2 Covl J f 
'1/2  1 


(s)ds  , fj(l/2) 


-19- 


Now,  Var  (?  j d/2))  was  computed  previously  and 


Var^f  « 


h (n)  1/2 


'ill  ,t' 

ii 

h(n)  , 


d s dt 


where  we  restrict  -i > 0 and  - *■  1 

n n 


Finally, 


u a 


2 Cov  J f (s)ds,  f.(l/2) 

\l/2  1 1 I 


2 • 0 


2 n 


fj  . 


v a.  K 


nh  (n)  j-0  •’ 


vi 

L 


/i-I  - 

J\ 

h(n) 


- K 


ds 


where 


1/2  , j = 0 or  n 
^1  , otherwise 


What  dors  this  converge  to? 


We  then  look  at 


nh(n) 


= L nb  <%r) 


ru  _JL_ 

* j / 2 »»(«>  \ h(n) 


which  converges  to  2a  J K (v)  dv 

-1 


2 nh(n)  Cov  1 J f (s)ds  , f (1/2) 
M/2  , 


►-20*"  J K“(v)  dv 


Thus 


(A  \ 2*2 
^2<U /”*  J*  K (v) 


5. 


Asymptotic  variance  when  X is  random. 


From  section  3.  3 , 


Thus  rQ2<u)  = rQj(u)  + I(u) 

A A . 

Var^rQ2(u)j  = Var^rQj(u)^  + Var^I(u)^  + 2Cov 
From  section  3.  1 , 

V«(rQl(»))  = ^ 02(Q(u>)  • / K2(«)  dt 

V"M  = 72  • iife)  Var  (m'  (X)  | X = Q(u)) 
4n 

where  m/  (x)  = r^-E[Y|X  = x] 

ox 

“d  2Cov(rQl(u).  I(u|)  a i • jjjL-  • C(u) 

It  then  follows  that 


nh(n)  Var(rQ2(u))-^a2(Q(u))  • J K2(t)dt  . 


To  compute  the  asymptotic  variance  of  rQ^(»)  , we  proceed 
again  by  steps  as  in  section  4.2  : 


(A  . J ^ A / ^ , 

rQ3(u)J  = Var^J  rQj  («)ds  + rQ^-y 


v«(rQl(i))  i 02(Q(i))  J K2(.)d. 

1 


/ U ^ V U U A A 

Var ( J rQ'(s)ds)  = J*  J Cov  (rQj'(s),  rQ/(t)Jdsdt 
1/2  1/2  1/2  V 7 

Cov^Q^s),  rQj'(t))  i 

-i-  | J J (uav  - ».)  K'(^)k'(^)  du  Q(»)  du  Q(v) 
nh  (n)  ^0  0 

t r o2(q(u|)  K'(i^i)  K'(S^)  du]  . 


UU  , A A UU 

J J*  Cov(rQ  '(s),  rQ '(t)J  d9  dt  = - J J C(s,  t)dsdt 
1/2  1/2  ' n 1 / 2 1/2 

♦ sr  / °2Mii>  (Kfef)  - ^w)}2<ix  • 


Finally, 


2 Cov(;  rQ/(a)do.  rQ  4))  = --r~r  02q4)  J K2(t)dt  + 


constant 


i 


6.  Preliminary  Conclusions 

A study  of  mean  integrated  squared  error  done  by  Melzer  (1978)  for 

A 

sample  sizes  n = 20,  50,  100  allows  us  to  conclude  that  rQ^f*)  does 

A 

not  improve  on  rQj(*)  . Also,  there  is  much  to  be  gained  by  normalizing 
the  estimators  so  that  the  weights  add  up  to  1 exactly.  This  has  no 
effect  on  our  asymptotic  results. 

A 

The  proposed  estimator  rQ^(*)  was  abandoned  after  a few  tries 
on  simulated  data  because  of  its  oscillating  behavior. 


-25- 


REFERENCES 


Clark,  R.  M.  (1977).  "Nonparametric  estimation  of  a smooth  regression 
function.  11  J.  Roy.  Statist.  Soc.  B,  39,  107-113. 

Melzer,  M.  (1978).  "Empirical  study  of  quantile  regression  estimators," 
to  be  published. 

Parzen,  E.  (1977).  "Nonparametric  statistical  data  science:  A unified 

approach  based  on  density  estimation  and  testing  for  'white  noise',  " 
Technical  Report  No.  47,  Statistical  Science  Division,  SUNY/ 
Buffalo. 

Yang,  S.  S.  (1977).  "Linear  function  of  concomitants  of  order  statistics,  " 
Technical  Report  No.  7,  Department  of  Mathematics,  Massachu- 
setts Institute  of  Technology. 


